Localization of classical waves in two dimensional random media: A comparison 
between the analytic theory and exact numerical simulation 
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The localization length for classical waves in two dimensional random media is calculated exactly, 
and is compared with the theoretical prediction from the previous analytic theory. Significant 
discrepancies are observed, ft is also shown that as the frequency varies, critical changes in the 
localization behavior can occur. Possible reasons for the discrepancies are discussed. 

PACS numbers: 43.20.+g, 42.25.Fx 



INTRODUCTION 

Propagation of waves through random medium has 
been and continues to be a subject of vivid research 
jl], ||, |[ |], fl ■ When propagating through media con- 
taining many scatterers, waves will be repeatedly scat- 
tered by each scattcrcr, forming a multiple scattering 
process Q. Multiple scattering is responsible for phe- 
nomena such as random laser ^ ^| , electronic transport 
in impure solids |lOfl , and photonic or acoustic bandgaps 
PI, |T^ . p"3[ . Under proper conditions, multiple scattering 
leads to the unusual phenomenon of wave localization, a 
concept introduced by Anderson [juj| . 

Such a localization phenomenon has been character- 
ized by two levels. One is the weak localization associated 
with the enhanced back scattering due to constructive 
interference from the reversed propagating paths. The 
second is termed as the strong localization (for brevity 
just termed as localization hereafter), in which signifi- 
cant inhibition of transmission surfaces, indicating that 
the energy mostly remains in a region of space in the 
neighborhood of the emission. The generally accepted 
wisdom on the connection between the weak localization 
and the strong localization is that enhanced backscatter- 
ing of a diffusive wave packet leads to an effective reduc- 
tion in the diffusion constant of the wave packet. When 
the influence of the increased backscattering is so over- 
whelming that the diffusion constant vanishes, the strong 
localization sets in. 

It is worthy to mention that tremendous efforts have 
been devoted to investigate the localization phenomenon 
for classical waves in random media over the past sev- 
eral years @ g || [gj [l| || [§, || g g|. How- 
ever, observation of classical wave localization is a diffi- 
cult task, partially because suitable systems are hard to 
find and partially because observation is often compli- 
cated by such effects as absorption and attenuation. 

One important quantity associated with the wave lo- 
calization is the localization length which is defined as 
the characteristic decay length of the transmitted inten- 
sity through the system from the source. Taking into 



account of the enhanced backscattering, the corrected 
(renormalized) diffusion constant may be obtained which 
then may be set to zero to obtain the localization length. 
This theoretical version of the localization length may 
be found in the literature for various dimensions for both 
the electronic waves and classical waves (e. g. [O, Q , and 
summarized in Ref. |p5[). 

A question we raise in this paper is about the validity 
of the existing theoretical formula for localization length, 
and subsequently the appropriateness or accuracy of the 
previous theory of wave localization. A main reason be- 
hind raising this question is that it would be helpful to 
experimentalists in designing experiments, provided that 
the theoretical formula for localization length is accurate 
enough in both qualitative and quantitative sense, and 
otherwise, it may mislead the experimentalists if the the- 
ory is not accurate enough. 

The main aim of this paper is to answer the questions 
raised above for classical waves. One possible way is to 
concentrate on an exactly computable system and find 
the exact characteristic decay length of the transmitted 
intensity with the sample size, and compare with the 
theoretical value of this particular system obtained from 
the existing theory. In this paper, we first outline the 
previous theory and show a numerically exact procedure 
for the calculation of the localization length for classi- 
cal waves in two dimensional media. Then we compare 
the two results, to gain information about the accuracy 
and appropriateness of the previous theory. In addition, 
we will also investigate the localization behavior as the 
frequency varies. These will be followed by a detailed 
discussion on the implications of the results. 



FORMALISM 
Existing theory 

Here we briefly review the existing theory for wave lo- 
calization. As wave propagates in random media, it ex- 
periences multiple random scattering, and as a result, the 
wave loses its propagating phase, leading to the gradual 
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decreases of the coherence of the wave in the absence of 
absorption, i. e. elastic scattering. Meanwhile, diffusive 
wave is built up as more and more scattering takes place. 
The traditional picture of localization is in fact a version 
of localization of diffusive waves. In other words, the con- 
ventional wisdom towards the localization mechanism is 
the absence of diffusion. The procedure to obtained the 
localization state can be briefly summarized as follows. 

The quantity, D^ B ' which is a measure of diffusion of 
classical waves is called the classical Boltzman diffusion 
constant and it may be derived under the coherent po- 
tential approximation , and is given as [^5| 



D (B) 



Vtl 

ddira 



(1) 



where Vt is the transport velocity, / is the mean free path 
and ddim is the dimensionality. 

As waves scattered along any two reversed paths in 
the backward direction interfere constructively, leading 
to the enhanced backscattering effect, which will then 
add corrections to the diffusion coefficient. In the two 
dimension case, such an enhanced backscattering effect 
is represented by a set of maximally crossed ladder di- 
agrams. An evaluation of these diagrams leads to an 
integration for which two cut-off limits have to be intro- 
duced to avoid the divergence. In this way, the correction 
to the diffusion constant for two dimensional system may 
be found to be IM 



SD ~ -\n(L M /l m ) 



(2) 



where Lm and l m are the two cut-off limits. It is then 
interpreted in the previous theory that the cut-off limit 
l m is a measure of the minimum scaling for the waves 
and is thought to be related to (for example) the mean 
free path, whereas Lm is a measure of the effective size 
of the sample. It is rather important to note that the 
correction in Eq. (||) is not only negative but diverges 
as Lm approach infinity. This is obviously unphysical, 
since the conductance or the corrected diffusion constant 
cannot be negative. To avoid the problem, it was con- 
jectured that Lm is in fact related to the localization 
range, or simply the localization length, in such way that 
when Lm is equal to the localization length denoted by 
£ say (i. e. Lm = the corrected diffusion coefficient 
becomes zero, i. e. the absence of diffusion: 



D R (0 = + 6D(0 = 0. 



(3) 



The localization length £ is solved from this equation. 
It is obvious that this equation always renders a solution 
and thus a localization length can always be found in two 
dimensions. Such an absence-of-diffusion mechanism is 
the core of the previous theory of localization in two di- 
mensions, and supports the assertion from a scaling anal- 
ysis that all waves are localized in two dimensions |2fl] . 



As will be shown later, there are significant discrep- 
ancies between the result from the absence-of-diffusion 
mechanism and the result from the exact numerical com- 
putation. These discrepancies would indicate either that 
the theory based on the diffusion mechanism is not accu- 
rate enough or simply that the theory has not yet fully 
captured the essence of the localization. We will come 
back to this later. Since the existing theory heavily de- 
pends on the diffusion picture, we may call it the diffusion 
based theory. 

Equation (|3|) leads to the following solution for the 
localization length in two dimensions 



theory = /exp(-Re[fc cff ]Z), 



(4) 



Under the effective medium theory, the mean free path, 
I may be obtained from the characteristic decay length 
of the coherent intensity and may be expressed as I = 
p-^j and the effective wave number, k c g may be writ- 
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ten as 



fcoff = k 




Pf(0) 



(5) 



where k is the wave vector of the incident wave prop- 
agating through the random scattering media with the 
scatterer numerical density p; /(0) represents the scat- 
tering function of a single scatterer in the forward direc- 
tion. The forward scattering function can be calculated 
by the standard series expansion method for a cylinder 
with radius a when a wave with wave vector k is incident 
normally on it. The forward scattering function is given 
as 



n=0 



(0) 



with 



C n — 



e n {J' n {ka)J n (k'a) - ^ j„(fc a ) J' n {k'a)) 



H^ ] \ka)J n {k'a) - ^H^\ka)J' n {^a) 



where e = 1 and e„>i = 2, g is the ratio of the den- 
sity of the scatterer to that of the media, h is the ratio 
of the wave speed in the sactterer to that in the media, 
J„ represents the Bessel function of order n, Hn repre- 
sents n-th order Hankel function of first kind and prime 
indicates the derivative. 



Exact numerical calculation 

Thus, on one hand the localization length for classi- 
cal wave in a system consisting of any kind of cylindrical 
scatterers placed randomly in any media may be evalu- 
ated by the use of equations from the previous theory. On 
the other hand, the localization length for classical waves 
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in two dimensional random media may also be calculated 
exactly by the use of the multiple scattering technique. 

Consider N straight cylinders located at f% with i = 
1,2, N to form a completely random array. An acoustic 
line source transmitting monochromatic waves is placed 
at f s . The scattered wave from each cylinder is a re- 
sponse to the total incident wave composed of the direct 
wave from the source and the multiply scattered waves 
from other cylinders. The final wave reaching a receiver 
located at r r is the sum of direct wave from the source 
and the scattered waves from all the cylinders. Such scat- 
tering problem can be solved exactly, following Twersky. 
While the details are in p?j| , the essential procedures 
are summarized below. The formulation presented be- 
low has been successfully applied to explain the recent 
experimental observations of acoustic crystals^, ^9j. 

The scattered wave from the j-th cylinder can be writ- 
ten as 



Ps(r,rj 



iirAiHW(k\r-r-- 



(7) 



where k is the wavenumber in the medium, is the 

n-th order Hankel function of first kind, and 4>r-fj is the 
azimuthal angle of the vector f—fj relative to the positive 
x axis. The total incident wave around the i-th cylinder 
(i = 1, 2, JV; i 7^ j) is summation of the direct incident 
wave from the source and the scattered waves from all 
other scatterers, can be expressed as 



B*J n (*|f-rJ|)e* 



(8) 



The coefficients A l n and B l n can be solved by expressing 
the scattered wave p s (r, fj ) , for each j ^ i, in terms of the 
modes with respect to the z-th scatterer by the addition 
theorem for Bessel function |3(J. Then the usual bound- 
ary conditions are matched at each individual scattering 
cylinder. This leads to 



B' 



N 

E c ^ 



with 



and 



(9) 



(10) 



CI' 



£ iirAjH^jm - r-D^-^-l , (11) 



and 



where T % n are the transfer matrices relating the acoustic 
properties of the scatterers and the surrounding medium, 
and have been given by Eq. (21) in Rcf. |2g| ]. 

The coefficients A l n and B 3 n can then be inverted from 
Eq. (||). Once the coefficients A l n are determined, the 
transmitted wave at any special point is given by 



P(r) =Po{r) 



N oo 

E E 

i—l n— — CK 



i7TAiH^(k\f-r-l\)e in ^ 



(13) 

where po is the field when no scatterers are present. 
The normalized transmission is defined as T = p/po and 
therefore, the acoustic intensity is represented as \T\ 2 . 

The averaged localization length is subsequently deter- 
mined by]3l[| 



< ln|T(i)| 2 >' 



(14) 



Bi =mriAl 



(12) 



where L is the sample size, and < • > denotes the ensem- 
ble average over the random distribution of the scatter- 
ers. Thus obtained localization length can be compared 
to that in Eq. ([|) obtained analytically from the previous 
theory. 



RESULTS COMPARISON AND DISCUSSION 
Results comparison 

The system we consider here consists of N identical air 
cylinders placed randomly in water medium. The reason 
behind considering the system as the air cylinders in wa- 
ter is that due to the large contrast of densities for air 
and water, and also a large contrast of sound speeds in 
air and water, the air cylinders act as strong scatters to 
the waves propagating in water media. The radius of 
each air cylinder is a. The fraction of area occupied by 
the cylinders per area is p. the average distance between 
nearest neighbors is therefore d = (^) 1 / 2 a, which is also 
the lattice constant for the corresponding regular lattice 
array. All the cylinders are placed completely randomly 
within a circle of radius L. A transmitting line source is 
located at the center and the receiver is located outside 
the scattering cloud. In the computation, the acoustic in- 
tensity is normalized in such a way that its value equals 
unity when there are no scatterers present. Thus the un- 
interesting geometrical spreading effect is naturally elim- 
inated. All the lengths are scaled by the parameter d, 
and the frequency is presented in terms of ka to make 
the computation non-dimensional. 

Fig. (1) shows the variation of < In \T\ 2 > as a func- 
tion of the system size, L for three different values of 
frequency (ka = 0.01, 0.013 and 0.016) with the filling 
fraction, (3 = 0.001. The < ■ > implies the configuration 
average of the total transmitted intensity. Number of 
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configurations considered here is 200. It is apparent from 
the Fig. (1) that the averaged total transmitted intensity 
decays exponentially with the dimcnsionless system size, 
L for all the ka values shown in the figure which in turn 
indicates the localization of wave at those frequencies. 
The localization length is nothing but the characteristic 
decay length and thus it may be obtained from the in- 
verse of the slope of the straight lines (in Fig. (1)) for 
those frequencies. Similarly one may obtain the localiza- 
tion length of classical waves at other frequencies from 
the exact numerical computation. 

Fig. (2) presents the dimensionless localization length 
(£) as a function of frequencies (ka) for (3 = 0.001. The 
dashed curve with circles represent the localization length 
obtained from the exact numerical computation using the 
multiple scattering technique while the solid curve rep- 
resents the localization length calculated from the the- 
oretical formula as prescribed earlier in Eq. (Q). From 
the dashed curve we observe that the localization length 
increases very rapidly as the frequency of the wave is 
decreased at the low frequency regime, the localization 
length also increases rapidly with the increase of fre- 
quency at the high frequency regime and there is a mini- 
mum of localization length at an intermediate frequency. 
For lower values of frequencies the waves propagating 
through the random media feel the media as homoge- 
neous over a large distance depending on the wavelength, 
for higher values of frequencies the waves are like very 
localized objects that can easily propagate between the 
scatterers and at an optimal frequency, the waves feel 
the strongest scattering in the random media. Thus, the 
larger values of localization lengths at lower and higher 
frequencies and the minimum value of localization length 
at an intermediate frequency are understood physically. 

However, we observe that critical change in the local- 
ization length takes place at lower and higher frequencies. 
On the other hand, the solid curve also shows similar be- 
havior. But, one important difference is that the mini- 
mum of the dashed curve occurs at a higher frequency 
when compared with that of the solid curve. In other 
wards, the strongest localization predicted from the pre- 
vious theoretical formula occurs at a frequency which 
is lower than that obtained from the exact numerical 
computation. And as a result, we observe that the two 
curves cross each other at some frequency, say (ka) cr . For 
ka < (ka) cr the localization length obtained from exact 
numerical calculation is larger than that obtained from 
the existing theoretical formula while for ka > (ka) cr , 
the theoretical value dominates over the exact value of 
localization length. [Note that both the results are the- 
oretical though, the results (values) obtained from the 
previous theoretical formula are termed as theoretical re- 
sults (values) and those obtained from exact numerical 
calculations are termed as exact results (values) and we 
may use this convention for our convenience.] In the low 
and high frequency regimes, the discrepancies between 



the theoretical and the exact values are significant. 

We have also done similar calculations for (3 = 0.01 
and compared the localization length obtained from the 
theoretical formula with that obtained from the exact 
numerical computation and is shown in Fig. (3). While 
the theoretical formula seems to capture the shape of the 
exact results (Fig. (2) and Fig. (3)), but there are sig- 
nificant discrepancies with regard to whether the waves 
are most localized, and where localization length starts 
to increase rapidly. 

Shown in Fig. (4) is the variation of localization length 
(£) as a function of the filling factor (3 for ka — 0.009 
which is less than (ka) cr in Fig. (2). The dashed curve 
with circle represents the exact localization length ob- 
tained from numerical calculation, and the solid curve is 
obtained from theoretical formula. The solid curve shows 
the monotonous decrease of the localization length as the 
filing factor increases but the exact result (dashed curve 
with circle) shows completely different behavior: the lo- 
calization length initially decreases and then increases 
monotonically, leaving a minimum at some value of f3. 
Thus, the wave with frequency ka = 0.009 becomes local- 
ized most strongly only for a specific value of the average 
density of scatterers. However, for ka = 0.009, the exact 
localization length is larger than the theoretical values 
for all filling factors. Thus here we observe that at the 
frequency, ka = 0.009 the theoretical results do not even 
capture the trend of the exact results. Similar behavior 
can be observed for a small range of frequencies around 
ka = 0.009. 

The variation of localization length (£) as a function 
of the filling factor (3 for a moderate frequency ka=0.015 
which is greater than (ka) cr in Fig. (2) is shown in 
Fig. (5). The dashed curve with circle represents the 
exact result and the solid curve represents the theoreti- 
cal result for localization length. The localization length 
decreases with filling factor for both the curves. However, 
there is a cross over between two curves at some value 
of 0, say j3 cr below which the theoretical value is higher 
than the exact value and above (3 cr the exact values are 
larger than the theoretical values for localization length. 
Thus at ka = 0.015, the exact and theoretical curves are 
similar in shape, the exact and theoretical values for lo- 
calization length matches at (3 = (3 cri otherwise there are 
significant quantitative discrepancies. All the observed 
discrepancies between the theoretical and exact results 
are probably due to the approximations involved in the 
previous theory. 

In order to study the fluctuation behavior of local- 
ization, we define an exponent, by analogy with one- 



dimensional cases) 32 



7 ^- I ln(|T| 2 ) 



(15) 



which is reciprocal to the localization length. Fig. 6(a) 
shows the variance of 7, i. e. var(7) =< 7 2 > — < 7 > 2 , 
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as a function of ka for [3 — 0.001. 

We interestingly observe that for a range of frequencies, 
the fluctuation is rather small indicating that the wave 
is strongly localized in this range of frequencies. How- 
ever, at low and high frequencies (below and above the 
range of frequencies where the fluctuation is very small) 
the fluctuation is high. This is an indication of some 
critical change in the localization behavior and it tends 
to suggest that waves become very rapidly and critically 
less localized at low and high frequencies. It may also 
be an implication of the localization-delocalization tran- 
sition, as discussed inj33|, [jij. These cannot be explained 
in the context of the current diffusion based theory. We 
also plot the variance of the exponent 7 versus the mean 
< 7 > in Fig. 6(b). Here we do not observe the linear de- 
pendence between the mean localization length and the 
variance, similar to ID cases|32, p5f 



Discussion 

All the deviations from the predictions from the ex- 
isting theory should possibly have a couple of important 
implications. The first is merely that the theory is not 
accurate enough, since it is necessarily based upon a per- 
turbation scheme given the complications involved in the 
problem. If this scenario is valid, we could believe that 
though inaccurate, the present theory still captures the 
fundamental nature of the localization phenomenon. It 
would then be expectable that the discrepancies will be 
reduced when more and more relevant multiple scatter- 
ing processes are included in the theoretical derivation. 
A next task would thus be to look for these scattering 
processes. 

The second implication would be simply that the exist- 
ing theory, which is based upon the absence-of-diffusion 
mechanism, has not yet fully encapsulated the essence 
of the phenomenon of wave localization. Beside the re- 
sults shown above, there are other evidences indicating 
that this line of reasoning should not be discarded too 
lightly. For example, it has been shown that in the local- 
ized state, a genuine phase-coherence behavior prevails 
as a unique feature of localized waves |B3]. This coherent 
behavior can be understood as follows. For quantum me- 
chanic or acoustic waves, the current can be written as 
J ~ Re[V>(— i)V^>], where ip stands for the wave function 
for quantum mechanical systems and for the pressure in 
acoustic systems. Writing the field as tp = \tp\e , the 
current becomes J ~ |^>| 2 V#. It is clear that when 9 
is constant at least by domains while \ip\ ^ 0, the flow 
stops, i. e. J = 0, and the wave or the energy is localized 
in space, i. e. \tp\ ^ 0. Obviously the constant phase 8 
indicates the appearance of a long range ordering in the 
system. All these have been demonstrated successfully 
not only for two dimensional media p3[ |34[ ], but for one 
and three dimensions as wellpll, |35|, pqj. It can be easily 



shown that the same consideration also holds for electro- 
magnetic waves ]37| . In the current diffusion based the- 
ory, the phase information is not attained. Therefore, at 
least the theory cannot account for the phase-coherence 
related localization phenomenon. Additionally, the the- 
ory could also be expected to fail when waves come to a 
complete stop before the diffusion becomes dominant. 

We would like to make a further note on the existing 
theory of wave localization. It has been widely believed 
in the literature that the theory has been tested exper- 
imentally to be very successful. Upon a close scrutiny, 
we find that the success is mainly based upon two types 
of experiments. One is the measurement of the en- 
hanced backscattering (e. g. Ref. JL8|). In a rigor- 
ous numerical simulation, it has been shown that the 
enhanced backscattering is not directly related to the 
localization |20| . Waves are not necessarily localized when 
a strong enhanced backscattering exits, and sometimes 
waves can be localized even when the enhanced backscat- 
tering is weak. Asides from few exceptions [p)[ 23 , the 
other type of experiments is based on observations of the 
exponential decay of waves as they propagate through 
disordered media, as summarized in 22 ; this setup sce- 
nario is in fact also the base for the previous scaling 
analysis f26f . The recent research p8|, however, has indi- 
cated that while this type of experiments can indeed help 
obtain information about whether a disordered medium 
will or not prevent waves from propagation through the 
medium, but this type of experiments is not sufficient 
to discern whether the medium really only has localized 
states; this discovery is consistent with the prediction of 
a recent scaling analysis J39|. Unwanted effects of non- 
localization origin can also contribute to the wanted ex- 
ponential decay, likely making data interpretation am- 
biguous. In other words, there is a crucial need to differ- 
entiate the situation that waves are blocked from trans- 
mission from the situation that waves can be actually 
localized in the medium. This is also essential to a scal- 
ing analysis 1 39 . It is worth noting that there was report 



of the observation of microwave localization in two di- 
mensions when a transmitting source is inside disordered 
media [Q. However, the diffusion based theory has not 
been verified against this experimental result. Instead, 
an exact numerical simulation based on the phase coher- 
ence picture is shown to agree well with the experimental 
observation 



SUMMARY 

The localization length for classical waves in two di- 
mensional random media is obtained from exact numeri- 
cal computation as well as from the previous theoretical 
formula and then the results are compared. Significant 
discrepancies between the two results is observed. Fur- 
thermore, results indicate some critical change in the lo- 
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calization behavior at low and high frequencies, and it 
seems to suggest that waves become very rapidly and 
critically less localized at these frequencies. These fea- 
tures are not readily seen in the previous theoretical pic- 
ture of wave localization. Therefore it is suggested that 
caution should be taken when applying the previous an- 
alytic theory to infer the localization length of classical 
waves. Possible implications of the present research are 
discussed. 
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Figure Captions 

Fig. 1 < ln|T| 2 > is plotted as a function of sample 
size (L) for ka = 0.010 (circles), ka — 0.013 (stars) 
and ka = 0.016 (diamonds) respectively. Here the 
filling factor (/3) is 0.001. 

Fig. 2 Localization length (£) is shown as a function of 
frequency (ka) for /3=0.001. The dashed curve with 
circles represent the exact values obtained numer- 
ically while the solid curve is obtained from theo- 
retical formula. 

Fig. 3 Same as the Fig. 2 except that (3 =0.01 here. 

Fig. 4 Localization length (£) is shown as a function of 
filling factor (j3) for fca=0.009. The dashed curve 
with circles represent the exact values obtained nu- 
merically while the solid curve is obtained from the- 
oretical formula. 

Fig. 5 Same as the Fig. 4 except that ka — 0.015 here. 

Fig. 6 (a) This shows the variance of the localization ex- 
ponent var(7) as a function of frequency (ka) with 
(3 = 0.001. (b) The variance var(7) is shown as a 
function of the mean < 7 > for [3 = 0.001. 



